Collapse of the vortex-lattice inductance and shear modulus at the melting transition 

in untwinned YBa 2 Cu30y. 
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The complex resistivity p(uf) of the vortex lattice in an untwinned crystal of 93-K YBa2Cu307 
has been measured at frequencies ui/2-r from 100 kHz to 20 MHz in a 2-Tesla field H || c, using a 
4-probe RF transmission technique that enables continuous measurements versus ui and tempera- 
ture T. As T is increased, the inductance C s (cj) = Im p{u))/uj increases steeply to a cusp at the 
melting temperature T m , and then undergoes a steep collapse consistent with vanishing of the shear 
modulus C66 ■ We discuss in detail the separation of the vortex-lattice inductance from the 'volume' 
inductance, and other skin-depth effects. To analyze the spectra, we consider a weakly disordered 
lattice with a low pin density. Close fits are obtained to pi(ui) over 2 decades in uj. Values of the 
pinning parameter k and shear modulus C66 obtained show that cm collapses by over 4 decades at 
T m , whereas k remains finite. 

PACS numbers: 74.60.Ge,74.72.Bk,72.15.Gd,64.70.Dv 



I. INTRODUCTION 

The investigation of the response of vortices in type 
II superconductors to an alternating driving current has 
had a long and fruitful history, beginning with the exper- 
iment of Gittlcman and Rosenblum (GR) on Pbln and 
Nb alloys Q. In cuprate superconductors, these inves- 
tigations have been performed on YBa2Cu307 (YBCO) 
over a broad range of frequencies extending from 100 kHz 
to a few THz. At frequencies u> below a few 100 MHz, 
the vortex response is obtained by directly measuring the 
sample's complex resistivity p(w). At microwave frequen- 
cies (10 to 100 GHz), cavity perturbation techniques have 
been the primary approach, although bolometric absorp- 
tion techniques have been more useful at the high mi- 
crowave end. 

The occurence of the vortex solid-to-liquid melting 
transition in the cuprates is one of the most interest- 
ing phenomena in the investigation of 'vortex matter' in 
superconductors [Q. The melting transition at T m has 
been investigated in some detail by high-resolution mea- 
surements of the magnetization |j, |J], heat capacity || 
and dc resistivity ^, 0). In any solid-to-liquid transition, 
the key quantity of interest is the shear modulus which 
vanishes in the liquid. The resistivity profile (p vs. T) in 
a crystal of Bi2Sr2CaCu20s in which alternating strong 
and weak pinning channels are created by ion irradiation 
was interpreted as a loss of shear strength at T m || . The 
effect of the vanishing shear modulus at T rn on p when 
the field H is tilted relative to the twin boundaries in 
twinned YBCO crystals was also demonstrated 

The standard approach to investigating how the shear 
modulus changes in a system undergoing a solid-to-liquid 
transition is by observing its dynamical response to an os- 
cillating driving force. The frequency dependence of p(u>) 
in the vortex state in YBCO has been investigated by 
several groups at frequencies in the 1 to 100 MHz range. 
Initial measurements using thin-film samples found that 



p(lo) = p\{ui) + jp2(u>) displays strong, featureless dis- 
persion at these frequencies JO], [l^, [l3). However, the 
response in high-purity single crystals is qualitatively dif- 
ferent. Wu, Ong, Gagnon and Taillefer (WOGT) (nj 
investigated the ac response in untwinned crystals of 
YBCO in the MHz range, and found that the inductance 
vs. H exhibits a step-wise change at the melting field. 
As a liquid cannot produce an inductive response, the 
step-wise change in the inductance is direct evidence for 
a collapse of the shear modulus cqq at T m . 

The 2-probe experiments of WOGT were performed at 
fixed T versus H . This approach cannot be extended to 
measurements in which T is continuously varied in fixed 
H because variations in the background signal (from T- 
dependent stray cable reactances) swamp the sample sig- 
nal. We report results obtained by the use of an RF 
4-probe technique that circumvents this obstacle. The 4- 
probe technique allows high-resolution impedance mea- 
surements versus each of the three quantities u>, H and T. 
We find that the collapse of the inductance is even more 
abrupt when T is scanned at fixed H. While we largely 
confirm the fixed-T results of WOGT, the swept-T ex- 
periments provide a more direct picture of the behavior 
of the vortex resistivity and inductance as T is increased 
above the melting line. With the broader perspective, we 
correct a previous inference regarding a 'remanent' shear 
modulus in the vortex liquid state. A detailed discussion 
of how the vortex-lattice inductance is separated from 
the 'volume' inductance is provided. Together with the 
measurements of WOGT, the present results provide a 
complete experimental picture of the linear response of 
the pinned lattice in a high-purity cuprate crystal to an 
ac driving current in the MHz frequency range. The two 
distinguishing features are the strong dispersion of p(u) 
observed even at low uj (just below the melting temper- 
ature T m ) and the striking abrupt collapse of the induc- 
tance at T m . 

To understand the spectra, we have adopted a model in 
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which the vortex lattice is weakly disordered by pinning 
to a low density of pins [lTj . A mean- field solution 
proposed by Ong and Wu (OW) Q is used to fit the 
measured resistivity spectra. By measuring the damping 
viscosity t](T) independently we reduce the number of 
adjustable parameters to just 2 numbers at each T. We 
show that the OW solution can achieve close fits to pi{to) 
extending over 2 decades in lj. The fits enable us to 
find both k and cqq (to within a multiplicative constant) 
at each T. The latter exhibits a remarkable 4-decade 
collapse to zero within 1 K of T m . 

In Sec. II a summary of the experimental approach 
is given. Section III discusses skin-depth effects and the 
'volume' inductance. In Sec. IV, we report measure- 
ments versus T and H at constant to. The central re- 
sults - the spectra of p\ and Pi)lj - are reported in Sec. 
Section VI explains the fits of the s pect ra to the 



weakly disordered lattice model, and Sec. |VII| discusses 
the temperature dependence of the shear modulus near 
T m . Appendix A describes the 4-probe technique, while 
Appendix B summarizes the impedance calculation for a 
sample with elliptical cross-section. Appendix C summa- 
rizes the linear response of the weakly disordered lattice 
to an ac current. 



of YBCO (of dimension a = 1.11 mm, b = 0.41 mm, and 
c = 55 /im), in which the dc resistivity displays a sharp 
transition (width of 0.2 K) at T c = 93.3 K. The RF resis- 
tivity was measured with the RF current density J || a, 
and H || c. Although all measurements were performed 
in a screened room, radio broadcast signals are picked 
up as weak resonance lines in the spectra. Parasitic line 
resonances, detected as weak resonances, are the hardest 
to eliminate. However, they can be minimized by empir- 
ically adjusting the grounding strap configuration. 

With the 4-probe method, we may vary independently 
each of the variables H (0-8 Tesla), lj (100 kHz to 20 
MHz), and T. To confirm that we are in the linear- 
response regime, we checked that curves of /5(w) mea- 
sured with incident power at -15.5 dBm and at 4.5 dBm 
are closely similar (except for a higher noise content 
in the -15.5 dBm curve for p^). We report all mea- 
surements in terms of the sample's complex resistivity 
p(u>) — Z s (uj)(bc/a') = pi(uj) + jp2(uj) 1 where a' = 0.68 
mm is the separation of the voltage contacts. In place 
of P2(w), it is preferable to discuss the (specific) induc- 
tance defined by C s = pijui (note that L s — L s (bc/a')). 
Dividing p2 by lj isolates the divergent behavior of C s as 
lj —* 0, thereby highlighting the dramatic collapse of C s 
at the lattice melting temperature. 



II. EXPERIMENTAL DETAILS 

As in WOGT [io| ], a weak, alternating current Je Jwt is 
applied within the Cu02 plane (J || a) in the presence of a 
field H || c. The average velocity response of the vortices 
is detected as a complex resistivity p(cj) — v(uj)B/ Jq. 

Because contact resistances on YBCO crystals are typ- 
ically a fraction of an Ohm, they add a T-dependent 
background that dominates the sample signal. In an ex- 
periment in which the field H is swept at constant T (as 
in WOGT), this large background may be subtracted, 
by reference to the zero-field data at each T. However, 
when T is swept at constant H, the contact background 
leads to substantial errors. The problem is exacerbated 
by the T dependence of the dielectric filling in the long 
coax cables. We minimized this by using special cables 
with low-density Polytetrafluoroethylene (teflon) dielec- 
tric. The higher propagation velocity (0.85c versus 0.66c) 
is also an advantage. 

The 4-probe design was developed to remove the back- 
ground altogether. The sample (impedance Z s = R s + 
juiLg) shunts the inner and outer conductors of both the 
incident and transmission coaxial cable (see Appendix 
A). If Z s <C Zq (the line impedance), the incident signal 
is strongly reflected, and the transmitted signal is highly 
sensitive to small changes in Z s . By phase-detecting the 
latter we may determine Z s . By inserting broad-band 
transformers in both the incident and transmission ca- 
bles, we can achieve true 4-probe measurement at fre- 
quencies up to 50 MHz (Appendix A describes the cali- 
bration procedure). 

The measurements were made on a detwinned crystal 



III. SKIN-DEPTH AND VOLUME 
INDUCTANCE 

In zero magnetic field, the impedance Z s {ui) of a super- 
conductor (in the form of a cylinder) may be calculated 
from the 2D Hehnholtz equation (V 2 + k 2 )A = 0, where 
k 2 = A~ 2 + 2j5~ 2 (here, A is the London penetration 
length and S the skin-depth determined by the quasipar- 
ticle conductivity). 

First, we consider the normal-state impedance (T > 
T c ) by letting A — > oo, and replacing 8 by the normal- 
state skin-depth S n = y / 2p n /ujpo (where p n is the 
normal-state resistivity and po the vacuum permeabil- 
ity). When lj is large enough that S n <C VA, J is con- 
fined to the skin-depth, and the effective R s (lj) increases 
(A is the cross-section area) . For a circular cross-section, 
Z s {lj) is given by §§ 



Z s (lj) 



Ljpo£S n Jq(x) 

2nr ^/2j Ji(x)' 



(1) 



where x = —jrQy/2j/6 n , J m {x) is the Bessel function of 
order m, and ro (£) is the radius (length) of the sample. 
The case of elliptical cross-section is treated in Appendix 
B. 

Equation |] shows that, at large lj, R s (lj) — > 
p£ I (2TTr 5 n ) , as expected. More important for our in- 
ductance discussion, the RF field energy at high u> is con- 
fined to the skin-depth. Therefore, the inductance L s (lj) 
decreases monotonically to zero as lj — > oo. As this in- 
ductance expresses the RF energy stored in the volume 
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of a straight cylinder, we refer to it as the 'volume' in- 
ductance L vo i(lu) to distinguish it from the vortex-lattice 
inductance which is the main focus of this work (and from 
stray, geometric inductances that arise from coil-coil flux 
linkage). From Eq. El we find that 



lim L vo i 



1^1 
2 4tt ' 



(2) 



For a high-purity Cu wire (rrj = 0.4 mm, 1= 3 mm) at 
300 K, L vo i(w) equals 120 pH at low u, but falls to a 
tenth of this value at 20 MHz. 

Below T c , the penetration of the RF fields (in zero H) 
is set by the London length, i.e. k — > A -1 . As the in- 
plane A in YBCO is very short (~ 0.2/xm), we may regard 
the stored RF field energy to be nominally zero. Hence 
the change in L s (in zero field) between T c and, say 60 
K, provides a direct measurement of L vo i. Applying Eq. 
^| to a cylinder with circular cross-section area equal to 
that in our crystal, we find that lim^^o L vo i — 34 pH 
(corresponding to C vo i = 1.2 x 10~ 13 Hem). This tiny 
inductance may be measured with about 2% resolution 
by our 4-probc method. 

In the resistive state (H 3> H c i), the RF penetration 
is determined by p\{uj) (i.e. k 2 — 2jSj 2 , where the flux- 
flow skin-depth 8f = \j2p\/uj^o) . As pi may be as large 
as 8 pilcm for H — 2 T, the volume inductance L vo i nearly 
recovers its full value given by Eq. [I]. However, in our 
measurements, the inductance increase is actually far in 
excess of L vo i. The excess inductance arises purely from 
the vortex lattice. 

Returning to the resistive component R s (u), we have 
converted it to p\ by the assumption that J is uniform 
over the cross-section (i.e. pi(co) = R s (uj)bc/a'). How- 
ever, if S is comparable to the half-thickness c/2, R s (u>) 
is significantly enhanced by the skin-depth effect, and pi 
will be overestimated. We now address this concern. As 
the aspect ratio b/c ~ 8.3, we cannot rely on Eq. [l] for 
quantitative guidance. In Appendix B, we have calcu- 
lated Z s (uj) for a normal-state sample with an elliptical 
cross-section using Mathieu functions. The variation of 
R s (us) and C vo i(u>) is shown in Fig. |l|) for a sample with 
this high ellipticity. We assumed the normal-state resis- 
tivity p n — 4/if2cm. 

As the results depend only on 5 n , values for R s (lu) and 
L vo i{io) at other values of p n are obtained by re-scaling 
the a;-axis. For example, at 85.9 K, pi measured in a field 
of 2 T with uj/2it — 3 MHz equals ~ 1/ificm. Re-scaling 
the frequency axis in Fig. [j] by a factor of 4, we see 
that at oj/2it — 12 MHz, R s (oj) and L vo i(u>) deviate by 
only 5% and 3%, respectively, from their low-frequency 
(i.e. uniform-distribution) values. At the lowest T (< 
84 K), however, the deviations may be up to 20% for 
uj/2tt >10 MHz, but these low T-high-w data will not 
play any significant role in our analysis. Throughout the 
Iow-oj region near the melting transition, the correction 
is less than 0.5 %, so that the uniform- J assumption is 
valid. 
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FIG. 1: Calculated impedance Z 3 (ui) = R a (u>) + juL s (uj) of 
a conducting cylinder of elliptical cross-section and resistivity 
p = 4/zficm. The inductance L s is entirely volume inductance 
L vo i. The parameters £o — 0.108 and h = 250pm are selected 
to match the sample cross-section. 



IV. FIXED-FREQUENCY MEASUREMENTS 

All the data reported here are taken either in zero field 
or in a fixed field of 2 T. In the resistive state, the induc- 
tance of YBCO varies strongly with u and T. We first 
display the variation of C s (u>,T, H) versus T measured 
at selected to in a 2-Tesla field as well as in zero field 
(main panel of Fig. ||) . The inset shows the correspond- 
ing curves for pi . 

In the main panel, we have plotted the specific induc- 
tance C s = P2/U with an arbitrary origin, because the 
total observed inductance L h s is the sum of the sam- 
ple contributions (L vo i and the vortex lattice term L v ) 
and an arbitrary background L\, g that arises from small 
uncertainties in setting 'null' in the phase setting of the 
lock-in amplifiers), viz. 



obs 



-L.M -I- Li; 



Jbg 



(3) 



We take the portion of the curve at 1 MHz-OT below 90 K 
as our reference, and identify sample inductance signals 
relative to this reference curve. 

We first discuss the 1 MHz curve in zero field. As T 
increases from 60 K, C s (u>i, T, 0) is nearly T-independent 
until, at T Cl it undergoes a step- increase to the normal- 
state value (wi/27r = 1 MHz). Clearly, the step corre- 
sponds to the collapse of the Meissner effect at T c , as dis- 
cussed above. The difference between C s (uii,T, 0) mea- 
sured at 100 and 60 K gives C vo i =1.0 x 10~ 13 Hem, 
rather close to the value 1.2 xl0~ 13 Hem estimated be- 
low Eq. H Hence we identify the step change in the 
curve C s (u>i , T, 0) with the maximum (specific) volume 
inductance C vo i of our sample. The curve C s (uii5,T, 0) 
measured at <jJi^/2tt = 15.4 MHz shows a closely similar 
step, as expected from Meissner flux-expulsion. 

When a field H = 2 T is applied, the sample induc- 
tance is dominated by that of the vortex lattice. We 
follow the curve at 1 MHz and 2 T, starting at 60 K. 
As T increases, C s (u>i,T, Hq) increases rapidly, rising to 
a cusp just below the melting temperature of the lattice 
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FIG. 2: (Main panel) Temperature dependence of the total 
sample inductance C s (lo,T,H) — p2(oj, T, H)/iv measured in 
a 2-Tesla field (H || c) at the frequencies indicated. Here, '15' 
means 15.4 MHz and H (in Tesla) is in parantheses. As T 
approaches T m from below, C s (measured at 1-4 MHz) rises 
to a sharp cusp before collapsing to the value C vo i observed in 
the normal state. For comparison, the zero-field traces taken 
at 1 and 15.4 MHz are also shown. The inductance scale has 
an arbitrary origin (because of the term L\, g in Eq. |^). The 
portion of the curve 1 (OT) below 90 K serves as our reference 
line. Inductance signals measured relative to this reference are 
identified as contributions from the sample. The step-change 
at T c in zero H reflects the Meissner expulsion of the stored 
RF field energy (reduction of volume inductance C vo i). The 
inset shows p\ measured at the same uo. 



T m ~ 88.5 K. Curves at higher frequencies show a similar 
but reduced peak. Above T m , C s collapses rapidly to a 
value very close to the normal state value. Because of the 
finite value of the RF vortex resistivity p\ (inset), expul- 
sion of the RF fields from the interior is now determined 
by the flux- flow skin-depth Sf = \j2pij 'pquj (Sf ^ A). 
At these temperatures, <5/ > c/2, so the RF fields are 
uniform in the vortex liquid state; there is no change in 
C s in crossing T c — 93.5 K. 

The sharp increase in C s as T approaches T m from 
below is one of our main results. The collapse of the in- 
ductance above T m is the constant- iJ-scan version of the 
collapse observed by WOGT in constant-T scans versus 
H. 

[As mentioned in Sec. I, the present swept-T scans 
provide a much clearer picture of how the inductance 
evolves into the normal-state value. In WOGT jl2), the 
finite value of C s above the melting field H m was inter- 
preted as evidence for a small, residual shear rigidity in 
the vortex liquid state (their Fig. lb). This inference 
is incorrect. As discussed here, the residual inductance 
above H m is simply the volume inductance C vo i . Within 
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FIG. 3: The high frequency resistivity pi(u>) measured at u> 
= 15.4 MHz versus applied field H || c at temperatures close 
to T c in untwinned YBCO (J || c). The minimum is caused by 
the 'peak effect' which appears close to the melting field H m . 
At each T, the initial slope of pi vs. H is used to determine 
the viscosity rj a (T). The T dependence of r/ a is plotted in 
the insert (r/;, obtained from a second sample is also plotted) . 
Near T c , r) a ~ (1 - T/Tc) 1 - 5 , while r) b ~ (1 - T/T c ). 



our resolution, the shear rigidity in the vortex liquid state 
is zero. This removes the need for the arbitrary 'back- 
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ground' liquid term p 2 introduced by WOTG 

At our highest frequencies, the lattice response is pre- 
dominantly resistive. As discussed below, the friction 
term in this limit dominates the lattice forces and pin- 
ning forces, so that pi is proportional to the reciprocal 
of the viscosity T]{T). Figure [| (main panel) displays 
the field dependence of pi(u>) measured at the frequency 
oj/2tt = 15.4 MHz at temperatures close to T c (the cur- 
rent is along the et-axis). At each T, we have fitted the 
initial slope to the (free) flux-flow expression Pf = B(j)Q jr\ 
in order to find rj(T) (plotted in the insert) [<f>o = h/2e 
is the flux quantum]. We measured f](T) in 2 untwinned 
crystals. In the geometry J || a, the uncertainty in deter- 
mining r\ a is large because the broad 'peak effect' restricts 
the range of H over which p\ vs. H is truly linear. More 
reliable results are obtained (for 77^) in the geometry J || b 
where the peak effect is much narrower. The two sets of 
viscosity r\ a and r\h are shown in Fig. || (insert). The 
viscosity data are used in the fits discussed in Sec. VI. 



V. FREQUENCY DEPENDENCE OF THE 
COMPLEX RESISTIVITY 

Figure displays the frequency dependence of 
pi(w,T, Hq) at temperatures from 80.6 to 88.9 K (H is 
fixed at 2 T). At temperatures T < T m ~ 88.5 K, the 
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FIG. 4: (a) The frequency dependence of pi(uj,T, Ho) at 
temperatures 80.6 to 88.9 K in a 2-Tesla field H || c. The 
curves for pi become increasingly u dependent as T m (~ 88.5 
K) is approached from below. In the vortex liquid state above 
T m , pi rapidly approaches a nominally u-independent value 
(curve at 88.9 K). (b) A subset of the data in upper panel 
plotted against logw to show the low-tj region. The solid 
lines are fits discussed in Sec. |vi]. 



value of pi in the limit lo — > is zero. However, above 
T m , the limiting value becomes finite. At the highest u, 
pi approaches the Bardeen-Stephen value for free flux- 
flow p f (inset of Fig. |J) . 

Below T m , pi is strongly dispersive. Close to T m , p\ 
displays a crossover from a sharply increasing region at 
low lo to a gradual region at high uj. The 'knee' feature 
separating the two regions rapidly moves to high frequen- 
cies as T falls below T m by a few K, eventually moving 
out of our frequency window. 

Figure || displays C s (u>,T, Ho) versus u> for tempera- 
tures below and above T m (with Hq = 2 T). As explained 
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FIG. 5: The frequency dependence of the sample inductance 
C s (uj,T, Ho) measured in a 2-T field Ho || c. £ 3 (referenced 
to the zero-field inductance curve) is the sum of C v and C vo i ■ 
The upper and lower panels are measurements taken below 
and above 87.95 K, respectively. As T — > T m from below (a), 
the vortex term C v increases sharply, especially at low ui. Just 
below T m (b), it collapses to zero, leaving the nominally flat 
spectrum that we identify as the volume inductance C vo i- 



following Eq. ^, we identify the sample inductance using 
the zero-field curve (below 90 K) as reference, viz. 

C,(u,T,H) = [L obs (u,T,H)- L obs (w, T, 0)] {be/ a') . 

C s rises to a sharp cusp as 
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As discussed in Sec. 
T —>■ T~, before undergoing a steep collapse above T m . 
The most pronounced dispersion occurs at low ui where 
C s shows divergent behavior in the limit to — > 0. 

Although the spectra of the inductance are not used in 
our fitting process in Sec. VI, we discuss various contri- 
butions to the observed inductance signal (which includes 
the volume inductance) for the sake of completeness. We 
may estimate C vo i(to, T, H) by appealing to the observed 
pi and the calculated curves in Fig. [|. Note that C vo i 
can never exceed 0.12 pHcm. At high frequencies, L v is 
negligible, so that the curve at 15.4 MHz (2T) in Fig. || 
may be taken as that for C vo i vs. T at large to (we ignore 
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the slight dip near T m ). Using the values of p\ measured 
at 15.4 MHz (Fig. [|), we find this interpretation is con- 
sistent with the calculations in Fig. El 

We may also use the observed values of p\ to estimate 
C vo i at low frequencies. At the temperatures of interest 
(85-90 K), pi is sufficiently large, even at 200 KHz, to 
satisfy the uniform-J condition, i.e. C vo i is close to its 
maximum value 0.12 pHcm. With increasing u, C vo i de- 
creases slowly to its high-w value as seen in Fig. |. For 
example, at 88.2 K (in Fig. ||b), our estimate of the vol- 
ume inductance spectrum is a curve that starts at 0.12 
pHcm at low ui and decreases slowly to ^0.06 pHcm at 
high lo. The divergent vortex term is easy to distinguish 
from this small, featureless background. 

This concludes the purely experimental part of the re- 
port. The spectra in Figs. Gj a nd |^, complemented by the 
measurements of WOTG jlfj, constitute a rather com- 
plete quantitative description of the linear-response of 
the vortex lattice in untwinned YBCO at RF frequencies. 
In the following sections, we interpret the measurements 
within a specific model, in order to extract the tempera- 
ture dependence of the vortex-lattice shear modulus. 



VI. FITS TO THE SPECTRA 

A. Discussion of the model 

In the model used by GR [BJ, the vortex equation of 
motion is 



r/u 



J X Zfi 



(4) 



where u is the vortex displacement from equilibrium, and 
77 and k are the damping viscosity and the Labusch pa- 
rameter, respectively, and B || z. Measurements of p in 
thin-film samples O, [T^] and in untwinned crystals pc| ] 
of YBCO are in strong disagreement with the Lorentzian 
response p = (B<f>o/r))[l — juj p /uj]~ l predicted by Eq. || 
(here uj p = n/rj). In particular, the dramatic dispersion 
observed when the lattice melting line is approached, in 
either constant-i? scans (WOGT) or in the constant-T 
scans reported here, cannot be reproduced by Eq. ^, 
even if one replaces n by an 'effective' parameter that 
depends on u> and B. 

In the context of the vortex glass jlj], Dorsey |lj| 
has derived scaling relations for the w-dependent conduc- 
tivity ct(u>) at the vortex glass-to-liquid transition, viz. 
a(u>) ~ (— iu>) a . In thin-film samples of YBCO, the mea- 
sured spectra do in fact display power laws in uj over 2-3 
decades (1-500 MHz) fl2"| . The situation in untwinned 
crystals is very different, however. The spectra reported 
here (and in WOGT |§) do not follow power laws at 
all. Further, the melting transition is abrupt and weakly 
first-order rather than the continuous transition of the 
vortex glass. Under the high-purity conditions found in 
untwinned YBCO crystals, the observed spectra are more 
consistent with that of a weakly disordered vortex lattice 
with random pins than a vortex glass (as we now show). 



In their theory of the depinned vortex lattice, Schmid 
and Hauger (l(| and Larkin and Ovchinikov |lj] repre- 
sented the rigidity of the moving periodic structure by 
the lattice force matrix Df m . As discussed below, the 
non-Debye spectrum of p(u>) and its rapid change vs. T 
in the vicinity of T m reflect a rapidly changing length- 
scale in the problem. To describe this, it is necessary 
to retain the q (wavevector) dependence of the Fourier 
transform D s (q) of D ; s m . (Clem and Coffey |l!| have 
shown that approximating D s (q) by a constant D merely 
reproduces a Debye-like spectrum.) With inclusion of the 
pinning forces, the equation of motion is [fl6| 



77m 



1,1' • Ul' 



J X Zfi 



(5) 



where the sites with pins are indexed by i (the pin- 
ning force is short-ranged). To match the strong disper- 
sion observed, OW assumed that the pin distribution is 
sparse Q (Rq >as, where R is the average pin sepa- 
ration, and ob — y^ifio/B the lattice spacing). Appendix 
C summarizes the solution of OW. 

In response to the driving current 3e JUJt , the magni- 
tude and phase of the vortex velocity (averaged over the 
sample) produces the ac voltage that determines p. At 
low u), the inter- vortex forces dominate the friction force 
nil, so the velocity response is predominantly inductive 
{loC v 3> pi). In the opposite limit of large w, nu is dom- 
inant. As C v — ► 0, pi approaches the free-flow value 
B<f>o/v- Thus C v is a measure of the average restrain- 
ing force on each vortex. It is largest if the restraining 
forces are dominant, but vanishes when viscous damping 
dominates. 

The main feature of the solution is that the lattice 
propagator G(R, w), which transmits information on the 
vortex displacement at site R to its surroundings, has an 
effective range Rg = as y cqq / nuj that varies as 1/ 
(Eq. ^6j). At low frequencies, Rg ^> as, which implies 
that the motion of any one vortex is strongly correlated 
with that of a huge number of neighbors, a subset of 
which arc pinned (the 2D correlation volume is ~ Rq)- 
Hence the velocity response is large and inductive. With 
increasing to, the average vortex is restrained by fewer 
and fewer pinned vortices as Rg shrinks. The in-phase 
(dissipative) component of the average velocity increases, 
while the inductance drops. When to exceeds the charac- 
teristic frequency 



^66 



47rc 66 a% 



E 2 ' 



(0) 



Rq falls below the average pin spacing Rq. The ma- 
jority of vortices now respond as if they are free, and 
pi rapidly approaches the Bardeen-Stephen free flux-flow 
value pf = B<fio/r] while C v decreases to zero. 

The crossover is observed as the knee in pi at T just 
below T m (see Fig. [|a). As T approaches T m from below, 
the rapid softening of Cqq causes this crossover to occur 
at progressively lower frequencies. Finally, in the liquid 
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state above T m , the vanishing of c^e implies that virtually 
all the vortices are decoupled from the pins: p\ equals pj 
at all ui in the RF range, while C v is zero. 



B. The fitting procedure 

As discussed in Appendix C, the OW solution has 4 
adjustable parameters at each T, viz. k(T), cqq(T), rj(T) 
and Rq. Empirically, the fits are not sensitive to the 
particular choice of Rq (as long as Rq ^> as)- Following 
WOGT, we set Rq = 7 a B for B = 1 T. 

As discussed in Sec. [V, the viscosity r](T) may be 
obtained from the initial slope of pi vs. H measured 
at 15.4 MHz (Fig. 0). Adopting the viscosity data in 
the insert of Fig. |3|7we then have only two adjustable 
parameters at each T, viz. k(T) and cqq(T) [equivalently, 
to p (T) and luqq(T)]. As we have a continuous spectrum 
measured over 2 decades in frequency to find 2 numbers, 
the fitting problem is evidently strongly over-constrained. 

The fitting procedure consists in selecting seed values 
for uj p and ujqs so that the quantity S(uj) (Eq. |2^) can be 
computed by summing the propagator G2_o(Ri,w) (Eq. 
p8| ) over an area much larger than Rq (the pins are ar- 
rayed in a triagular superlattice of spacing Rq, and the 
sum is typically extended to 12 impurity 'shells' from 
the origin). With S(uj), we may compute the complex 
resistivity p(u>), viz. 



p(uj) = 



<1 



1 - 



1 



(7) 



The calculated curve is then compared with the mea- 
sured spectra of p\(u), and the fit parameters are 
changed in an iterative manner until the fit is optimized 
(the inductance data are not used in the optimization). 
For spectra near T m , the high-frequency curvature is 
most sensitive to lo p whereas the curvature at low oj is 
largely determined by the choice of ujqq, as expected on 
physical grounds (the value of 7/ mainly fixes the over- 
all scale of pi). As shown in Fig. |Jb, close fits to the 
data are achieved over the 2-decade range of u>. At low 
l>, the calculated p\ is capable of matching the curvature 
of the measurements to a surprising degree (with only 2 
adjustable parameters). The close agreement achieved in 
Fig. ^b offers encouraging support for the validity of the 
model. 



C. Calculated inductance 

With the parameters ui p and W66 optimized at each T, 
we may calculate the vortex-lattice inductance C v . As 
discussed above, in the limit lu — > 0, the extended range 
of the propagator G2d(R,u) implies that each vortex 
is restrained by a great many pins. Hence the velocity 
response is purely inductive. The calculated C v displays 
the strong divergence observed in the RF measurements. 




12 3 4 

Frequency (MHz) 



FIG. 6: Comparison of calculated p2{uj)/ui (solid lines) with 
the measured vortex inductance £ v (u,T, Ho) (open symbols) 
over the frequency range 100 kHz to 20 MHz in a field Ho = 
2 T. pi(yj)juj is calculated from Eq. (fusing the values of lo v 
and W66 derived from the fits to pi in Fig. Wb. The estimated 
volume inductance has been subtracted from the total sample 
inductance (see text). Above 87 K, C v strongly diverges as 
w — » 0, while at 83.9 K, the divergence is not observed. 



In Fig. ||, we compare at 4 selected T the calculated C v 
(solid lines) with the observed C v (open symbols). 

The sample inductance is the sum of the vortex-lattice 
term C v and the volume inductance C vo i, viz. 

T, H) = C v {u;, T, H) + £ vo1 (lu, T, H). (8) 

The volume term cannot be measured directly at low w 
because it is much weaker than the vortex term. How- 
ever, at frequencies above ~ 3-10 MHz (depending on 
T), it accounts for nearly all of the sample inductance 
(see Sec. 



IV). Because of this uncertainty, we used the 



spectra of p\ to optimize the fit parameters. 

The contribution of the volume inductance has to be 
estimated to isolate the vortex- lattice term (Eq. ||). In 
Sec. IV, this was carried out using the observed p\ to- 



gether with calculations for a sample with elliptical sec- 
tion. The estimated C vo i has been subtracted to isolate 
the vortex term C v in the plots displayed in Fig. ^. Al- 
though the uncertainties in the estimated C vo i are largest 
at low o>, they have the least impact on the comparison 
because C s diverges steeply in this limit. The compar- 
ison (inset) shows good agreement between calculation 
and measurement above 87 K (no further refinements of 
the fit parameters were made in comparing with C v ). 
As T approaches T m ~ 88.5 K, the low-w divergence in 
the observed C v progresses to lower frequencies (see the 



curves at 87.9, 88.2 and 88.4 K in inset). The calculated 
curves also match this progression. At low T, however 
(83.9 K curve), the agreement is not as good because the 
uniform-J assumption used to extract p is increasingly 
suspect. 



VII. FIT RESULTS 
A. Collapse of the shear modulus 

The values of k(T) and cee(T), obtained from u> p (T) 
and uj e6 (T), respectively, are shown in panels (a) and 
(b) of Fig. |t]. In each case, we show two sets of values 
depending on whether rj a or rjh has been used in the fits. 
We interpret the decrease of k as T — > T c (Panel (a)) 
as reflecting the decrease in condensation energy. We 
note that k remains finite at the melting temperature T m . 
Hence our fit strongly argues against models in which the 
transition at T m is interpreted in terms of a vanishing k. 

The critical parameter is the shear modulus cqq which 
undergoes a steep collapse by more than 3 decades within 
an interval of 1 K. The collapse occurs over the back- 
ground (gradual) decrease of cqq associated with the 
power-law decrease of the condensate density as T — > T c . 
Because of the rapid collapse in cqg, the steep increase 
in C s at temperatures below T rn is interrupted. Within 
a 1-K interval, C s drops to its normal-state value (back- 
ground curve in Fig. |7)d). 

Each of the discrete data points for cqq is derived from 
a full spectrum for p\. To follow the collapse in more 
detail, we have also adopted a different procedure, using 
the 1 MHz curve for C s and p\ measured continuously 
vs. T, together with interpolated values of rj a and rjb. 
This allows us to determine cqq as a continuous curve 
very close to T m (the fitting procedure becomes unstable 
at lower T). We show these continuous curves as solid 
and broken lines in Fig. 0b. At low T, the continuous- 
fit values are slightly larger than the spectra-based fits 
which are more accurate, but close to T m , they match 
very well. The continuous curve allows us to monitor the 
decrease of cqq over 4 decades. 

From Fig. ^, we conclude that the transition at T m 
is associated with a rapid collapse of the shear modulus. 
However, the Labusch parameter n remains finite. This 
provides firm evidence that the transition involves the 
collapse of the shear modulus, rather than the vanishing 
of the Labusch parameter. 
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FIG. 7: (a) Values of the Labusch parameter re obtained from 
the fits, using either n a (open circles) or rjb (solid squares). 
The sharp peak (right scale) is the behavior of the sample 
inductance C s = P2/UI measured at 1 MHz in a 2-T field [the 
zero of C s is arbitrary (Sec. |Tv|)]. (b) The normalized shear 
modulus C66 (log scale) versus T obtained from the fits, using 
either n a (open circles) or r]b (solid squares) in Fig. ^. Near 
T m , C66 decreases by 4 orders of magnitude within an interval 
of 0.5 K. The solid line represent values of C66 obtained from 
fits to p(u) with lu fixed at 1 MHz, and H=2 T, using r/b- The 
dashed line is obtained with lu fixed at 2 MHz, but using 77a. 



B. Summary 

In summary, we have performed measurements of p{ui) 
of the vortex lattice in an untwinned YBCO crystal to 
investigate the collapse of the lattice inductance at the 
melting temperature, using a high-resolution 4-probe RF 
technique. At all T investigated, the resistivity pi(cj) 
is strongly dispersive. Over the frequency range 100 



kHz - 20 MHz, pi increases from zero to the free-flow 
value Bipo/r), with a cross-over frequency scale (toee) that 
rapidly decreases towards zero as T approaches the melt- 
ing temperature T m from below. In the solid phase, the 
inductance pijui displays a steep divergence as 10 — > 0. 
However, at T m , this low- frequency divergence collapses 
to give an ^-independent inductance in the vortex-liquid 
state above T m . The observed spectra are qualitatively 
different from those in thin-film YBCO, and incompati- 
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ble with the predictions of vortex glass theory. To extract 
the vortex-lattice shear modulus from these spectra, we 
have used the mean-field solution of the Schmid-Hauger- 
Larkin-Ovchinikov model to fit the complex resistivity 
spectra. From the 2-parameter fit at each temperature, 
the normalized shear modulus cqq (T) displays a remark- 
able 4-decade collapse towards zero at T m . Hence the 
observed collapse of the inductance is driven by the van- 
ishing of the shear modulus rather than the vanishing of 
the pinning strength k. 



Appendix A: Experimental details 

We describe first the 2-probe method used by WOGT 
p0| . The sample Z 2 bridges the inner and outer con- 
ductors of both the incident and transmission coax ca- 
bles of line impedance Zq (Fig. ||a). The load resistors 
Z\ and Z[ are 100-O thin-film resistors. This configura- 
tion maximizes the sensitivity of the transmitted wave to 
small changes in Z 2 when \Z 2 \ <C Zq. The incident wave 
is transmitted with transmission coefficient Tt-, and re- 
flected with reflection coefficient Tr, where 



2Z0Z2 



Z 2 (2Z + Z, + Z[) + (Z + Z X )(Z + Z[) 
Zq + Z[ Zi — z a 

1 t -f- . 

Zq + Z\ Z\ + Zq 



7v (9) 



(10) 



Phase-sensitive detection of the transmitted signal al- 
lows Z 2 to be determined (Fig. ||a). The incident power 
from a synthesizer (Hewlett Packard 3336C) is typically 
fixed at +4.5 dBm, corresponding to a current less than 
5 mA at the sample. The transmitted signal is phase- 
detected by a high-frequency lock-in Extender/Enhancer 
(Palo Alto Research PAR100), which down-converts sig- 
nals in the range 1-25 MHz to 1-50 kHz, suitable for a 
standard lock- in amplifier (Stanford Research SR830). 

To modify the technique for 4-probe measurements 
(Fig. ||a), we inserted an RF transformer Tl of band- 
width 10 kHz - 150 MHz (Mini- Circuits T1-6-KK81) at 
the output of the synthesizer. The transmitted signal 
passes through a second transformer T2 before entering 
the PAR100. In Fig. |8|b, the contact impedances are 
represented by the impedances Z S i and Z' si . The trans- 
mission and reflection coefficients arc now 



26 



(i + £t)(i + ^)-£ 2 2 ' 

(fr -!)(! + &)- £2 



(l + 6r)(l + &) 



(11) 

(12) 



where £ T = {Z x + Z sl + Z 2 + Z s3 )/Z , & = + Kx + 
Z 2 + Z' s3 )/Z , and &=Z 2 /Z . 

While the 4-probe method provides much higher reso- 
lution, a drawback is that the inserted transformers intro- 
duce strong reflections. The largest reflection comes from 
the upstream transformer Tl. A wave reflected from the 
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FIG. 8: (a) Schematic diagram of the measurement circuit. 
The sample shunts the inner and outer conductors, so that 
the transmitted signal is highly sensitive to slight changes 
in the sample impedance Z2. Transformers Tl and T2 are 
inserted to decouple grounds in the 4-probe technique. The 
lower figure (b) shows the labelling of the contact impedances 
Z 3 i, Z' si , which may greatly exceed Z2 in amplitude. 



sample is reflected again at Tl and adds an cj-dependent 
contribution to the original transmitted wave. 

The effect of the multiply reflected wave on the total 
transmitted signal at the detector is expressed as 



j^obs 
1 T 



Vo, 



e) k ^T T T T [1 + e? 2kh r R T R 



(13) 



where T T and T R are, respectively, the transmission and 
reflection coefficients of Tl. The path- lengths from Tl 
to the detector, and from Tl to the sample, are called 
Itot and t\, respectively. Higher-order reflection contri- 
butions (notably from T2) are indicated by • • • . In ad- 
dition, we are now also sensitive to slight w-dependent 
deviations from unity of the PARlOO's transfer function 

To compensate for these two background contribu- 
tions, we used a reference combination Z ri Zi <r , Z[ r tai- 
lored to have a reflection coefficient T r R nearly equal to 
that of the sample (hereafter, rfj). As \Z S \ -C 1, the sam- 
ple's Ffj is very close to 0.33, with Zi and Z[ — 100 f2. 
The best choice for the reference combination is Z r = 10 
O, and Zi jT = Z[ = 90.9 f2, which has an w-independent 
r^j of 0.3344, and a transmission coefficient T T = 0.044. 

With the sample in place, the signal presented at the 
input of the SR 830 lock-in is given by V a = g(u))V ou t (Z s ), 
where V ou t is given by Eq. |l3L If the sample (and load 
resistors) are replaced by the reference combination, we 
have V r = g{u>)V ou t(Z r ). Dividing these two equations 
removes g(uj). This leaves the multiple reflection factor 
(quantity in [•••] in Eq. |T^ ). However, as T S R and T r R 
are nearly identical by design, the leading term in [■ • • ] is 
the same (the other terms are down by a factor T T « 1). 
Hence we have 



v r r T 1 



0{T S 7 



(14) 
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where and Tj, are given by Eq. [fl] with Z 2 — Z s and 
Z r , respectively. Our procedure is to measure and store 
the curve of V r vs. u> at each temperature of interest. 
Then V s is measured versus to. Since is known, 
(and Z s ) may be calculated from V s and V r using Eq. 
[L4j . We checked the reliability of the procedure using a 
second reference consisting of just a single 50-J7 resistor 
connecting the inner conductors of the two coax cables 
(this also has T r R = 0.333, but a much larger Tip = 0.667). 



Appendix B: Impedance of a conductor with 
elliptical cross-section 

The impedance Z s (u>) of a conductor with elliptical 
cross-section is obtained by solving the 2D Helmholtz 
equation (V 2 + k 2 )A = in elliptical coordinates (£,rj), 
where k 2 = 2j5~ 2 . With A(£, 77) = if>(£)(()(ri)z, the 
Helmholtz equation separates into the two Mathieu equa- 
tions, viz. [El] 



d 2 ^ 

i. 

dij 



-j^- - [a{q) - 2qcosh2£]i; = 0, 



_, , [a(q) - 2q cos 2iq](j) = 0, 



(15) 
(16) 



where q = —k 2 h 2 /4 (h is the foci spacing). With the 
boundary condition A(£o,rj) = As (i.e. a constant on 
the surface S), the vector potential may be expanded in 
terms of the Matthieu functions as 



Appendix C: The solution of Ong and Wu 

We summarize the equations in Ref. |l8| . The dis- 
placement ui of a vortex at site 1 is described by the 
equation 



77U1 



^Df m • u m + k Ui4i = J x z0 o . (19) 



pins 1 



Only a subset of vortices (at sites indexed by i) are af- 
fected by the pins. The remaining vortices are 'free', 
apart from being restrained by the lattice forces repre- 
sented by Dj m (s =T,L indexes the polarization mode). 
The displacement at an arbitrary site 1 is related to the 
displacement of the pinned vortices by 



ui(w) 



J 4*o 



N 

e s (q)e s (q)e 



pins 1 

jq-(Rl-Ri) 



[D s {q)+iur,} 



(20) 



where e s (q) is the unit vector for the mode s. 

At the mean-field level, the displacement of the vortex 
at a pin-site Ui mp is 



-1 J^o 
iWpM = y - — 

1 



■£(u>y 



(21) 



where all the interaction effects are in the term 



K [G(0, W ) + ^G(R ; ,c)], (22) 



(17) 



n=0 



Ce 2n is expressed by the expansion Ce 2n (^,q) = 

J2T=o ^2r cosh(2r£) [the same expansion holds for ce 2n , 
but with cos(2r?7) in place of cosh(2r£)]. 

The Mathieu functions arc not tabulated because the 
separation parameter a(q) and the coefficients A^J 1 ^ are 
functions of q. For a given q, a(q) satisfies the continued 
fraction a = -(<f /2)/[l - (a/4) - (g 2 /64)/[l -•••)]] 
We truncated the continued fraction at the 12 th order. 
The 12 roots {a n (q)} generate the full matrix A^ 2 ™' with 
(n,r) = 1,- • -12. 

Finally, with the ac voltage V(ui) = —jujAs and the ac 
current I = — k 2 //jQ d£d?7A(£, 77), the impedance may 
be expressed as 



Z s (u>) 



(2tt) 2 



E 



(Ag n) ) 2 Ce' 2n (€o,q) 
Ce 2n (£ , q) 



-1 -1 



(18) 



where Ce' 2n — dCe 2n /d^, and C 2n =/ Q ' 1 ce 2n (ri, q)drj. In 
Fig. [I], we display the ui dependence of R s {lo) and C{u) 
calculated with parameters £0 — 0.108 and h = 250/jm, 
matching the sample cross-section, and p\ = 4^Slcm. 



and Up = n/r] is the pinning frequency. The lattice prop- 
agator G(R, to) is defined as Jl6[, |l8f 



G(R,w) 



1 v (y-6 s (q)) 2 e^ R 
[£> s (q)+jwr/] 



(23) 



In Eq. |2^, the first term represents the uniform q=0 
mode of the transmitted elastic forces. The remaining 
terms are elastic forces transmitted to the pinned vor- 
tex at the origin arising from its own displacement 
(G(0,u>)), or the displacement of other pinned vortices 
G(Ri,cj) (the prime indicates that the i = term is 
left out). 

Neglecting the longitudinal modes, the transverse force 
matrix is expressed as flq] 



D T (q) 



[ceeq 



2l 2 

Cii q z \a B . 



(24) 



where q is the in-plane component of q, viz. q = (g, g 2 ), 
and C66 and C44 are the shear and tilt moduli, respec- 
tively. 

In the 2D limit, the propagator simplifies to the form 
(after angular-averaging within the ab plane and neglect- 
ing the tilt modulus C44) 



r m ^ [ Q ,qMqR) 
kG 2D Ri,w =9 / d q , 2 , ■ 2i 

J a [q + JP \ 



(25) 
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where Q = y/Aw/ag, Jo(x) is the zeroth order Bessel 
function, R = |Rj_|, with R = (Rj_, R z ), and g the cou- 
pling constant k/Attcqq. The frequency dependence ap- 
pears only in the characteristic wavevector 



with [3 = 0.52 and 7 = 0.85. Equations |27] and £8] were 
given incorrectly in OW [Q (cp/ their Eqs. 18 and 19). 



P 



rjoj 



(26) 

Acknowledgments 



which serves as a cut-off of the logarithmic divergence in 
the 2D case. Setting R = in Eq. r2B, we have 



kG 2D (0,w) 




j arctan (- 



As the propagator G2d(R,w) is used repeatedly in the 
sum in Eq. |22j, it is convenient to adopt the approxima- 
tion (accurate to a few percent) 
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